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ABSTRACT 

Ab initio quantum chemical techniques have been used 
to investigate covalent ly- bonded and hydrogen-bonded 
species that may be important intermediates in the 
reaction of hydroxyl and hydroperoxyl radicals. 

Stable structures of both types were identified. 

Basis sets of polarized double zeta quality and large 
scale configuration interaction wave functions have 
been utilized. Based upon electronic energies, the 
covalently bonded HOOOH spec ies is found to be 26.4 
kcal/mol more stable than the OH and HO^ radicals. 

Similarly, the hydrogen bonded HO HO.., species is 

found to have an electronic energy 4.7 kcal/mol below 
that of the component radicals, after correction is 
made for the basis set superposition error. The 
hydrogen bonded form is found to be planar, to 
possess one relatively "normal” hydrogen bond, and to 
have lowest energy ^A ' and ^A' states that are 
essentially degenerate. The 1 A" and ^A” excited 
states produced by rotation of the unpaired OH 
electron into the molecular plane are found to be 
very slightly bound. 


PACS numbers: 31.20.Tz, 82.30.Cf, 31.20.Ej, 82. 40. We 
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I. INTRODUCTION 

Radical-radical gas phase reactions have received 
considerable experimental * r . J theoretical attention in recent 
years. In large part this has been motivated by their importance 
in atmospheric and combustion chemistry. In the past, it often 
proved attractive to view those reactions that resulted in simple 
atom transfer (metathesis) as proceeding via direct abstraction 
mechanisms. However, as detailed kinetics studies have become 
available, it is now clear that a significant number of these 
"simple" reactions actually display complicated kinetic 
behavior, 1 such as pressure dependence of their rate constants 
and negative activation energies. This is frequently 
rationalized in terms of mechanisms that proceed through 
stabilizable intermediate species. 

The hydroxyl and hydroperoxyl radicals play important roles 

2 

in the chemistry of both the lower and upper atmospheres. The 

OH radical is an important initiator of atmospheric hydrocarbon 

oxidation. Both the generation and destruction of ozone in the 

troposphere are partially accounted for by chemistry involving 

the HO species, and the coupling of HO and HO^. to the NO and 
x 1-0 2 x 

Cl^ cycles is very important in determining the loss rate of 
stratospher ic ozone. The concentrations of the OH and HO^ 
radicals are controlled to a significant degree by 
radical-radical recombination reactions. 

The reaction between the hydroxyl and hydroperoxyl radicals 
HO ♦ H0 2 --> H 2 0 ♦ 0 2 


( 1 ) 
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has received a great deal of attention because of its importance 

in atmospheric HO^ chemistry.^ There are now several direct 
4 — 6 

measurements of the bimolecular rate constant in the low 

pressure regime that are in agreement, giving values in the range 
- 1 1 3-1 

(6. 4-7. 5) x 10 cm s at total pressures of a few torr. At 

pressures of about one atmosphere, there is general agreement 

7-1 1 

among several experiments that k 1 is in the range 

(1.0-1. 2) x 10 1 ^ cm^s 1 . This apparent pressure dependence is 

included in the value of k^ currently recommended by the NASA 

1 2 

panel for use in atmospheric modeling. Regardless of its 
detailed pressure dependence, the rate of Reaction (1 ) is 
unusually fast. The super f ic ial 1 y similar self-recombination 

reactions of OH and HO^ are both much slower, with rate 

12 -123-1 

constants in the range (1-2) x 10 cm s Finally, a recent 

study 10 shows that k^ has a negative temperature dependence that 

can be described by an E /R of -416 in the range 250-420 K. 

a 

The data discussed above has led several authors to suggest 
that Reaction (1 ) may proceed via formation of an intermediate 
species. It is also possible that multiple mechanisms, both 
direct and indirect, are important. Thus far, there have been no 
direct observations of an intermediate, and such an experiment 
promises to be difficult. The goal of the present project is to 
apply ab initio quantum chemical methods to the study cf possible 
intermediates for Reaction (1 >. The structures of these 
intermediate species will be elucidated and their stabilities 
will be estimated. These stabilities are calculated as the 
energy of the formation reaction 
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OH ♦ H0 2 --> H 2 0 3 , ( 2 > 

in which H 2 0^ represents any proposed intermediate species. 

In determining which of the ^0^ intermediates are possible 

it is necessary to consider the electronic states of the 

2 

reactants and products. The reactant ground states, 1T( OH > and 
2 

A"(H0 2 ), car correlate with either singlet or triplet 

3 _ _ 1 

intermediates. The products have I ( 0_ ) and A„(H^0> ground 

g 2 12 

states, correlating only with triplet intermedi ates . However, 

the reaction is exothermic, and it is necessary to consider 

product excited states as the initially formed species. Using 

1 4 

standard heats of formation, the enthalpy change for Reaction 

( 1 ) at 298 K is calculated to be -69.6 kcal/mol. This energy is 

1 , 1 ♦ 

sufficient to open channels leading to the A and I excited 

g g 

states of 0 2> These excited state products can be correlated 
with a singlet i ntermedi ate , making possible both singlet and 
triplet spin-allowed mechanisms. Restrictions on the spatial 
symmetry of the intermediate exist only if the reaction is 
constrained to proceed via a path with high symmetry. 

15 IS 

It should also be noted that there are two studies ’ in 
the literature that have questioned whether the pressure 
dependence of Reaction (1 ) is real and, therefore, whether there 
is any need to invoke a mechanism with an intermediate species. 
In fact, if the isotope study of Reference 16 could be verified 
by other experiments, it would appear to place a serious 
limitation on the acceptable nature of any proposed intermediate 
While an intermediate in Reaction (1 > has not been studied 
directly and no species of the formula H 2 0 3 has been observed in 
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the gas phase, such a species has been generated and studied in 

17-19 


condensed phase experiments 


The vibrational spectra of 


both H^O-j and have been observed and bands have been 


assigned to the oxygen framework modes, assuming an open-chain 


HOOOH structure. 


1 8c 


The ultraviolet spectrum of H^O^ reported 


to be similar to, but more intense than, that of H^O^,. 


1 7b 


Early theoretical studies of the HOOOH molecule included SCF 


calculations using both minimal and extended basis sets with 


20-22 23 

partial geometry optimization. A recent SCF study 


employing minimal and split valence basis sets treated the HOOOH 


isomer and three possible hydrogen-bonded forms of H^O^. These 


workers reported that only the HOOOH isomer is stable relative to 


the HO and HO^ fragments. This is in disagreement with the 


results of the present study and will be further discussed in the 
conclusion section of this paper. The only previous study to 

24 

include a treatment of electron correlation is that of Cremer. 
This work employed basis sets ranging from minimal to polarized 
double zeta and accounted for the effects of electron correlation 
with second order many-body perturbation theory. In addition to 
the minimum-energy geometry of HOOOH, an extensive 

characterization of the potential surface as a function of the OH 
rotational angles was reported. 

The stability of the HOOOH intermediate relative to the OH 


and H0 2 radicals, i.e. the energy of Reaction (2), has not been 


accurately calculated from ab initio methods or measured 


25 


directly. From group additivity principles Benson has 


estimated the heat of formation of the H^O^ open chain form to be 


v-J 
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1 4 

-17.7 kcal/mol. Combined with tabulated heats of formation for 
HO (9.3 kcal/mol) and HO^ <2.5 kcal/mol), this yields an enthalpy 
change at 298 K for Reaction (2) of -29.5 kcal/mol. 


page 8 


II. COMPUTATIONAL METHOD 
A. SCF method and basis sets 

The molecular geometries were optimized by SCF gradient 

calculations that were carried out using the method and computer 

26 

program due to Komornicki . These programs use the unrestricted 

Hartree Fock method for open shell systems. As a measure of the 

UHF spin contamination , we note that the expectation value of the 
2 

S operator was no greater than 0.76 for the doublet cases and 

2.02 for the triplet case. The basis set used in these SCF 

«« 


gradient calculations was the standard 6-31C polarized split 
27 

valence set, except that the polarization functions were chosen 

to have exponents of 0.85 on oxygen and 1.0 on hydrogen. 

In all phases of this .study except the gradient 

calculations, a polarized double-zeta ( DZP ) quality Cartesian 

Gaussian basis set was used. In particular, all energies reported 

in this paper were calculated using the DZP basis set, which 

20 29 

consisted of Dunning's C4s2p|2sD contraction of Huzinaga's 

(9s5p|4s> set. A scale factor of £= 1.2 was chosen for the 

hydrogen s-functions, and polarization was provided by sets of d 

functions on oxygen (o=0.85) and p functions on hydrogen (cx=1.0>. 

Except for the gradient calculations, the open-shell SCF 

solutions were obtained in the spin restricted form (RHF), and 

the molecular orbitals were constrained within D_. subgroups to 

2h 

reflect the symmetry of the nuclear framework. 

B. Configuration interaction calculations 

°0 

Conventional Cl techniques with configuration selection 
and extrapolation were used. The master configuration list in 
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each case consisted of all single and double excitations relative 

to either a single- or multi-configuration reference set such 

that the oxygen Is-core orbitals remained doubly occupied. 

Generally, all single excitations were retained and the 

perturbation theory technique of Gershgorn and Shavitt^ 1 was used 

to select the most important double excitations. The selection 

procedure accepted or rejected a particular space configuration, 

with its entire set of spin couplings, based upon the sum of the 

A k contributions for that set. Beginning with the least 

important one, configurations were rejected in sequence until the 

sum of the energy contributions for the rejected conf igurations 

32 

was equal to T, the cumulative selection threshold. Cl 
energies E( x ) were obtained at three values of T , fit to a 
straight line, and extrapolated^ to zero threshold. The 
resulting energies are denoted as "CI(SD>" or "extrap CI(SD)". 

To obtain improved Cl convergence, approximate natural 
orbitals (ANO's) were used thruugnout. At each conformation, 
preliminary Cl calculations in the SCF orbital basis were 
performed to define the ANO's. Conf igurat ion selection was used 
to reject all conf igurations with estimated energy contributions 

-4 34 

less than 10 hartree. After freezing those orbitals that 

35 

were occupied in the SCF wave function, the natural orbitals of 
these small Cl expansions (1000-2500 terms) were obtained. The 
ANO's consisted of the occupied SCF orbitals plus the set of 
"lightly occupied" natural orbitals. The three ANO's with 
smallest occupation numbers, generally equal io zero, were core 
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correlation orbitals and were discarded, since the Is core 

orbitals remained fully occupied throughout this work. 

The most important effects of quadruple excitations were 

36 

accounted for in an approximate way using Davidson's correction 
f ormula , 

AE = ( 1 - C 2 )A E , , ( 3 ) 

q O d 

where AE . is the correlation energy due to the entire set of 
d 

double excitations, AE is the analogous quantity due to 

q 

quadruple excitations, and C is the coefficient of L.. in the 

O bL r 

"all doubles" Cl, In this work, the quantity AE, was estimated 

d 

to be equal to the extrapolated "singles plus doubles" 

correlation energy, and C was taken from the Cl calculation at 

o 

the smallest of the selection thresholds. Where a 

2 

multi-conf iguration reference state was used, C q was taken as 

the sum of the squares of the coefficients corresponding to the 

37 

reference configurations. These Cl energies that include the 
Davidson correction are denoted as "CI(SDQ)". 

C. Computer programs 

In the course of this work several conputer programs were 
used. These include the SCF gradient program of A. Komornicki, 
the MOLECULE integral program due to J. Almlof , the NASA Langley 
SCF-CI package due largely to C. W. Bauschlicher and B. H. 
Lengsfield, and the MELD package developed by E. R. Davidson and 


coworkers . 
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III. RESULTS AND DISCUSSION 
A. HO and HO^ radicals 

A primary goal of this study is the calculation of relative 
stabilities of possible intermediates for Reaction (1 >. In 
practice this is accomplished by calculating the energy change 
for Reaction (2), in which H^O^ represents the particular 
intermediate being considered. Care must be taken to ensure a 
well balanced treatment of the reactants and products so that the 
calculated energy difference will be an accurate estimate of the 
reaction energy. To this end, the energies of the OH and HO^ 
radicals have been obtained using the same basis set and Cl 

procedures that are used for the H^O^ complexes. 

38 39 

The literature ’ contains geometry optimization studies 
of both HO and HO^ using basis sets and Cl procedures that are 
very similar to the ones used here. Since the theoretical 
geometries for these species have been found to be in good 
agreement with experiment, reoptimization was considered 
unnecessary. Except as noted otherwise, the experimental HO^ and 
HO geometries are used in this work. For the OH radical the bond 


^ U 0 ] 

length is 0.971 A, and for HO^ the parameters are 0.977 A 

( Rqh > • 1*334 A(Rqq), and 104.2° (bond angle). The reference 

( SCF ) configurations for the ground states are 

2 


OH : 


, 2 _ 2 _ 2 , 3 
1 0 2o 3 a 1 - 


n 


22222222 ' 
H0 2 : la* 2a* 3a* 4a’ 5 a' 6a* la" 7a’ 2a” 

Equivalence restrictions were not imposed on the 


2 A" 


and 


( 4 ) 
( 5 ) 


orbitals of OH. 
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The energies obtained at three levels of theory for the OH 

and H0-, radicals and their sums are entered on the first three 

42 

lines of Table I. As is well known , truncated Cl expansions 

are generally net size extensive. In the present case, this is 

manifest in the fact that the sum of Cl energies for OH and HCu, 

at a particular excitation level will not be equal to the 

corresponding energy obtained for them as a supermolecule . On 

line four of Table I results are given for the supermolecule 

treatment of the H0_ and OH radicals at a separation of 500 a . 

d o 

This calculation was carried out for the (in C g symmetry) 

coupling of configurations (4) and (5). Comparison of the third 

and fourth rows of Table I shows that the SCF energies are size 

extensive as expected, but at the extrapolated C I ( SD ) level there 

is an error of 25.2 mh . Since the presence of a nonneg 1 i g i bl e 

size extensivity error means that higher excitations are 

important, it is expected that the Davidson correction for 

quadruple excitations will reduce this error. As the last column 

of the table shows, the error is reduced considerably at the 

C I ( SDQ ) level, but is still significant, with a value of 11.1 mh. 

It is somewhat difficult to estimate a priori the accuracy 

of the Davidson correction for a particular problem. One 

approach is to repeat the calculation using mu 1 1 i -con f i gurat i on 

reference states (HREF), thus including many of the most 

important quadruple excitations explicitly. This should ave a 

smaller correction to be made with the mul t i -ref erence analogue 
37 

of Eq <3> . In this spirit, MREF calculations were carried out 

using four-configuration reference states for both OH and H0 o and 
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a three-configuration reference state for the supermolecule 
calculation. The suppl ement ary reference configurations were 
those that had been found to be the most important in the 
s i ng le- re f er enc e < SREF > calculations. The ANO’s were not 
redefined at the MREf level. If the Davidson correction is 
reasonably accurate, the SREF and MREF results should compare 
favorably at the C I ( SDQ > level. In fact, the MREF value of E< OH > 
♦ E(HO_,> was only 0.5 mh below the SREF result, and the MRFF 
supermolecule energy was 1.3 mh lower than its SREF counterpart . 
These results suggest that the Davidson correction yields 
reasonably wel 1 -converged C I ( SDQ ) energies at the SRFF level. 

In the limit of very large reference sets, the size 
extensivity error should vanish. However, it was reduced only 
slightly from 11.1 to 10.3 mh by the present procedure, 
indicating that the supermolecule approach is still required at 
this limited MREF level. For this particular electronic state of 
the separated fragments the MREF and SREF results agree closely, 
and the MRFF description does not seem to be required for 
"chemical accuracy". That, of course, would not necessarily be 
true for other choices of geometries and electronic states. I t 
is concluded that C I ( SD ) with the Davidson correction yields 
results for this system that are converged to approximately i ml., 
and that the superinolecule treatment of the separated f ragmen’s 
should be used to account for the size extensivity error inherin’ 
in the Cl approach. 


i 


I 
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B. Covalently-bonded HOOOH intermediate 

1. Equilibrium geometry 

One possible intermediate in Reaction (1 1 results from 
formation of a covalent bond between the HO^, and OH radicals and 
is represented here by the formula HOOOH. The equilibrium 
conformation of this covalent molecule has not been determined 

exper imenta 1 1 y , but has been predicted by the theoretical study 

24 ** 

of Cremer . This previous work employed a 6-31G basis set 

that is of similar quality to the DZP basis used here and treated 

electron correlation via manv-bodv perturbation theory. It was 

expected a priori that the difference between Cremer 's ^2^3 

structure and one determined by the present Cl methods would be 

small. Since these structures would both correspond to the same 

minimum in the potential surface, the energy difference between 

them should also be quite small. This consideration could have 

led us to directly employ Cremer 's structure in these 

calculations . However, in order to make sure that there was no 

24 

inconsistency between his approach and our own, a partial 
geometry optimization of the H p 0^ molecule was carried out. 

The SCF gradient procedure converged to a conformation 
similar to the expected one, with the hydrogens above and below 
the 000 plane in an arrangement with symmetry. The electronic 
ground state is of ^ symmetry, and the SCF conf igurat ion is 
given by 


£.<* 


1b 3a 2b 4a 3b 


ba 


4b 6a 


7a“ 6b 


('A 


) C 6 > 
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A sketch is shown below: 


H 

f 

O^-o^o 

H 

The bond lengths and angles at the SCF level differed from the 

24 

previous results in a way that is consistent with the use of an 

uncorrelated wave function, with the following differences being 

noted: R ( - . 026 A), R (-0.067 A), 000 ang le( ♦ 0 . 9° ) , H00 
UH UU 

ang le( * 2 . 6° ) , and dihedral ang le( *8 . 0° ) . 

Using the DZP basis set and Cl wave functions, these 

molecular parameters were optimized sequentially. For this 

purpose we employed a cumulative Cl selection threshold of 75 mh 

without any extrapolation. After one complete cycle of 

optimization, the discrepancies with Cremer 's results had been 

reduced to: R QK <-0.004 A), R OQ (-0.03 A), 000 angle (+0.9°), H00 

ang 1 e( * 2 . 2° ) , and dihedral ang le< ♦5 . 8° ) . At this point it was 

clear that the geometry eventually obtained with the Cl method 

24 

would be rather close to the published result; certainly, the 

two methods appear to be consistent. From the magnitudes of the 

energy changes during the optimization cycle, it was estimated 

that further refinement would result in an energy drop of no more 

than 1.5 mh. Therefore, further optimization cycles were not 

carried out, and Cremer 's geometry was used for the remaining 

calculations: R_„(0.972 A), R (1.439 A), 000 ang le( 1 06 . 3° ) , H00 
UH UU 

ang le( 1 00 . 2° ) , and dihedral angl e( 78 . 1 ° ) . 
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2. Wave function and molecular properties 

The electronic structure and bonding in HOOOH has previously 

24 

been discussed by Cremer. We have examined the natural 

orbitals and find that they generally support his analysis. The 

rationalization of the OH bond orientation at the equilibrium 

geometry was given both in terms of dipole-dipole repulsion and 

43 

in terms of the electronic anomeric effect. . This latter 

argument is usually presented in a way that requires the lone 

pair on the terminal oxygen atom to be coplanar with the 

nonadjacent 0-0 bond, such that the lone pair may interact 

favorably with the 0-0 antibonding orbital. Examination of the 

contour plots for the lone pair orbitals on the terminal oxygen 

atoms showed that they are oriented within 5-10 degrees of the 

000 plane, making that explanation plausible. 

44 

The Mull iken populations for the largest Cl wave function 
show a considerable transfer of charge from the hydrogen atoms to 
the terminal oxygens, with only a little charge accumulating on 
the central oxygen. The electronic populations are 8.31 on 
terminal oxygen atoms, 8.07 on central oxygen, and 0.65 on each 
hydr gen atom. The electric dipole moment is constrained by 
symmetry to lie along the axis and points from the central 
oxygen atom (negative end) toward the plane of the hydrogen 
atoms. The OH bond moments are situated such that there is 
considerable cancel lat ion , giving a relatively small net moment 
of 1.26 Debye at the Cl level. 


3. Energies 
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In Table I are presented HOOOH total energies and the 
binding energies for dissociation into the OH and HC^ radicals. 
The geometry of the HOOOH molecule was that given above, and ("he 
extrapolated Cl energy was obtained using selection thresholds of 
100, 75, and 50 mh. At the 50 mh threshold the Cl expansion 
consisted of approximately 15 600 terms, and yielded a 
variational energy (before extrapolation) of -226.0752 hartree. 
The binding energy of HOOOH is calculated at each level of theory 
as the difference between Lines 4 and 5 of Table I. From the 
C I ( SDQ ) data this binding energy is ♦42.1 mh, or +26.4 kcal/mol. 
At the CI(SD) level this quantity is only +32.4 mh and at the 7CF 
level it is -16.5 mh (unbound). Thus, some treatment of electron 
correlation is essential in order to obtain the correct sign for 
the binding energy, and a careful treatment of its contribution 
is necessary for an accurate value. 

The CI(SDQ) electronic binding energies (Table I) are the 
"best” purely theoretical results presented here, in the sense 
that no experimental data or jad hoc corrections have been 
included. However, to facilitate comparison with experiment, it 
is desirable to include the zero point energy correction and the 
effects of temperature, as well as to estimate the inherent error 
due to basis set limitations. This resulting series of "derived” 
energies is discussed in subsequent paragraphs and presented in 
Table II. 

The thermodynamic energy change of Reaction (2) at 0 K, AE 

o 

can be calculated as the negative of the electronic binding 
energy plus a contribution from the zero point vibrational 
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energy. For this purpose all vibrational modes were treated as 

harmonic, and experimental values of the vibrational frequencies 

were used where available. Otherwise, the necessary frequencies 

were estimated by analogy to hydrogen peroxide. The vibrational 

frequencies of OH (3735 cm 1 and HO^ (3436, 1392, and 1098 
- 1 45 

cm ) were available in the literature. Of the nine 

vibrational modes in H^O^, the stretching and bending modes of 

1 8c 

the oxygen framework were taken from the matrix study (855, 

755, and 500 cm 1 ), while the remaining modes were estimated from 

46> 

the s P ec ^ ruiT1 * The OH stretches were both taken as equal 

to 3607 cm 1 , the corresponding value in H^O.^. Similarly, the 

-1 

0-0-H bending frequencies were set equal to 1330 cm , the 

average of the values, and the hindered rotations were both 

taken equal to the 317 cm 1 mode of hydrogen peroxide. Using 

these wave numbers, the difference in zero point energies of the 

product and reactants in Reaction (2) is 1478 cm 1 , or 4.23 

kcal/mol. Inclusion of this correction results in a calculated 

energy ( AE ) of Reaction (2) at 0 K of -22.2 kcal/mol. 
o 

The calculated energy of Reaction (2) includes the effects 
of higher order Cl excitations in an approximate way and 
represents a reasonably thorough treatment of the correlation 
energy contribution. The largest source of error remaining is 
probably due to the limitations of the basis set. In general, 
basis sets of medium quality such as the present one ( DZP ) 
describe the atomic situation more adequately than the molecular 
one. Consequently, bond dissociation energies calculated with 
such basis sets are usually too small if the computational 


page 19 


procedure is otherwise well balanced. When, as in the present 

case, the dissociation process goes to fragment molecules rather 

than atoms, the error may be somewhat smaller, but still 

significant. It was not feasible to repeat these calculations 

with a larger basis set. However, the literature contains 

several studies with basis sets of similar quality in which 

calculated binding energies are compared to experimental 

dissociation energies. These results have been employed to 

obtain an approximate correction to the binding energy. 

In the recent literature there are several reports based on 

MBPT or the closely related CCD approaches that attempt to 

calculate dissociation energies with basis sets of roughly the 

same quality as the present one. Since the treatment of electron 

correlation in the present Cl approach and that of the MBPT or 

CCD methods are quite comparable, it is expected that the 

conclusions about basis set deficiencies drawn from those studies 

47 

will be applicable to this work. Bartlett and coworkers have 

studied several systems with DZP basis and found the calculated 

47b 

binding energies to be too small by 1-6%. Adams et al. state 
that for this same quality basis set, one can expect calculated 

dissociation energies for single bonds to be 1-3 kcal/mol too 

** 48 

small. Using the 6-31 1C basis set, Pople and coworkers find 

dissociation and atomization energies to be in the range 4-15% 

too small if directly calculated for cases such as the present 

one, in which the number of unpaired electrons is not conserved. 

** 

Using both 6-31G and DZP basis sets, dissociation energies of 

49 

several first row diatomics have been calculated accurately to 
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within 5-152. The general experience is that calculated 
dissociation energies are almost always too small and are 
somewhat worse for the cases that involve breaking of multiple 
bonds or dissociation into atoms rather than into fragment 
molecules . 

The error estimates given above have broad ranges and 

exceptions can easily be found. However, it is clear that the 

calculated energy change for Reaction (2) is almost certainly too 

small in magnitude because of basis set limitations. We do not 

anticipate an exceptionally large error because Reaction (2) 

represents the breaking of a single bond. However, the error 

will probably be somewhat larger than in the very best cases, 

47b 

which seem to be the breaking of a bond involving hydrogen. 

Accordingly, we have elected to make an ad hoc correction of 62, 

or 1.3 kcal/mol. This gives a corrected value for the energy of 

Reaction (2), AE* of -23.5 kcal/mol. In this discussion the 

o 

prime on the thermodynamic quantities indicates that this 
correction has been included. 

For comparison purposes it is useful to derive an estimate 
of the room temperature enthalpy of Reaction (2) from our 
results. To convert AE ' to AE' it is noted that Reaction 

O c!70 

(2) represents the loss of three translational and two rotational 
degrees of freedom. Treating them classically, they contribute 
-2.5 RT (-1.48 kcal/mol at 2^8 K) to AE*. The vibrations are 
treated as harmonic oscillators and the internal energy above the 
zero point is given by the usual formula for an oscillator of 


frequency v 
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E . . = Nk 0/C expC P/T )-1 U 

vib 

where 0 is hv/k. Using the product and reactant wave numbers 
cnosen for the zero point energy calculations, the vibrational 
contribution to AE'^g is ca ^ culated to be ♦0.73 kcal/mol. 
Including the vibrational <♦0.73 kcal/mol) and 

rotat ional-translat ional (-1.48 kcal/mol) contributions at 298 K 
yields an energy of reaction AE' OQ of -24.3 kcal/mol. To 

t 70 

convert this to an enthalpy change, the work term (APV) is 
estimated as -RT (-0.59 kcal/mol), giving AH’ QQ = -24.8 
kcal/mol . 

25 298 

As stated in the introduction, Benson's estimate of AH^. 
of H^O^ implies a room temperature change for Reaction (2) of 
-29.5 kcal/mol. The disagreement between that value and the 
present one is a little less than 5 kcal/mol, a number that 
possesses at least marginal significance. Many of the 
uncertainties in the present work, such as the Cl extrapolation 
error, accuracy of the quadruple excitation correction, and 
errors in choice of unknown vibrational frequencies are 
essentially random and should combine with some cancellation. 

The estimate made of the error in reaction energy due to basis 
set limitations could easily be too small by as much as 2 
kcal/mol, but is unlikely to be more than 0.5 kcal/mol too large. 
If it were 2 kcal/mol larger, the present estimate and Benson’s 
would be in reasonable agreement. Rather than saying that these 
two numbers are in disagreement, it is more realistic to say that 
our results suggest that the correct enthaply change of Reaction 
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(2) may lie between the present value and that derived f rom 
25 

Benson's result. 


1 


w 
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C. Hydrogen-bonded Complex 

1. Minimum-energy conformation 

In addition to the covalently bonded HOCOH species, the 
various hydrogen bonded forms should also be considered as 
possible intermedi ates for Reaction (1). In this section we 
report the results of a potential energy surface study of the 
hydrogen-bonded form of H^O^ and a theoretical estimate of the 
energy of reaction for its formation 

H0 2 ♦ OH > HO H00 (7) 

Since hydrogen bonding does not involve large scale 
reorganization of the fragments’ electronic distributions, it is 
expected that the SCF description will be reasonably correct. 
Accordingly, the geometry optimization was carried out at only 
the SCF level, but the final energies were determined by both the 
SCF and Cl procedures. 

The minimum-energy conformation for the complex was located 
with the open shell ( UHF ) variant of the SCF gradient program. 

The electronic conf iguration used represents the triplet spin 
pairing of the OH and HO^ radicals. The minimum energy 
conformation was found to be planar and to possess a single 
hydrogen bond with the H0 2 radical serving as proton donor . The 
electronic wave function possessed C g symmetry and is represented 

by the conf iguration 

2 2 2 2 2 2 2 

la’ 2a’ 3a* 4a’ 5a’ 6a' 7a ' 

2 2 2 2 2 11 3 

8a’ 9a’ la" 10a’ 11a’ 2a" 3a" ' A’> (8) 

The unpaired electrons on both radicals are in oxygen p-type 

orbitals situated perpendicular to the molecular plane (a"). The 
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3a” MO, which represents the HO^ unpaired electron, is mostly 
localized on the terminal oxygen atom. A drawing of the complex 
is given in Fig. 1 along with the structural parameters. 

The structure in Fig. 1 strongly resembles unperturbed OH 
and HO^ fragments held together by one fairly "normal", albeit 
bent, hydrogen bond. Since the structure of the complex was 
determined at the UHF level , it should be compared to structures 
of the free radicals that were also calculated at that level. 
Therefore, UHF gradient calculations were carried out for the HO 
and HO^ radicals. When compared with these structures, it is 
seen that the bond lengths and angles within the fragments have 
undergone the following distortions upon complex formation: 

OH: R qh ( -0.002 X ) 

H° 2 : R oh (*0.003 X), Rqq( -0.002 X), H00 AngleC ♦ 1 . 4° ) . 

The effects of including electron correlation should be almost 
identical for the parameters of the isolated radicals and their 
complexed counterparts . Therefore, the distortions upon complex 
formation that are given above would be expected to carry over 
almost unchanged to geometries determined at the Cl level. 

2. Energy of the complex 

The binding energy of the hydrogen-bonded complex has been 
determined at the various levels of theory and the results are- 
presented in Table III. Cl calculations were carried out for the 
separate OH and H0 2 radicals and for their supermolecule using 
the UHF con f ormat ions . The size extensivity error as calculated 
from Lines 4-5 of Table III differs by only 0.2 mh from that 
found using experimental conf ormations (Table I). From Table III 
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it is seen that the SCF calculations predict the complex to be 
stable with a binding energy of 7.4 mh or 4.6 kcal/mol. At the 
C I ( SDQ > level this quantity increases to 8.5 mh or 5.3 kcal/mol. 

/. 1 though the correlation energy contribution to binding is only 
about 1-2 mh, it is significant because the binding energy itself 
is quite small. The use of either the supermolecule approach or 
an explicit size extensivity correction is essential in this 
situation, because otherwise, the Cl results would predict the 
complex to be unbound. 

In theoretical treatments of weak hydrogen bonds, such as in 
the case of the water dimer, some correction is often made for 
the basis set superposition error ( BSSE ) . This error results 
from basis set deficiencies and manifests itself in binding 
energies that are too large. If the basis set is inadequate, an 
improvement in the description of each f raement may be obtained 
,pon complex formation by making use of the orbitals of the other 
fragment. This improvement is an artifact of the basis set and 

has nothing to do with the real binding energy of the complex. 

50 

The counterpoise correction, a common method of estimating this 
effect, is determined by recalculating the fragment energies in 
the presence of the basis functions of the other fragment. Any 
improvement in energy is considered an estimate of the BSSE. 

This error can be very large when small basis sets are used. 

In the present case the BSSE at the SCF level was found to 
be 0.7 mh for OH and 0.3 mh for H0 Q , giving a total error of 
1.0 mh (0.6 kcal/mol). At the C I ( SD ) and CI(SDQ) levels the 
error was reduced by about 502. This reduction was unexpected 
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and may be due to fortuitous cancellation by the Cl extrapolation 

error, which is also on the order of 1 mh. Even allowing for 

such a cancellation, this result can be taken as evidence that 

the Cl BSSE in this case is not very much larger than that at the 

SCF level. Accordingly, we assume a BSSE correction of *0.6 

kcal/mol, giving corrected binding energies of 4.0 and 4.7 

kcal/mol, respectively , at the SCF and CI(SDQ) levels. 

The relatively small BSSE calculated in this case, 

especially at the Cl level, should not be taken as a measure of 

the completeness of this modest basis set. The magnitude of 

these errors depends strongly on the particular nature of the 

basis set being used and on the molecular conformation being 

considered. For example, the basis set used by Clementi and 
51 

Habitz in their water dimer calculation war certainly more 

complete than the present one, and their BSSE was smaller at the 

SCF level. However, at the Cl level the apparent error grew 

significantly. The reader is referred to Ref. 51 for a careful 

discussion of the Cl superposition error and the difficulty and 

52 

ambiguities involved in correcting for it. Also, Bauschlicher 

has pointed out that a large BSSE can arise with DZP quality 

basis sets when they are augmented with bond functions. 

3. Wave function and molecular properties 

In hydrogen bond formation the electrostatic interaction is 

53 

generally very important, and between polar neutral species the 
dipole-dipole interaction is its leading component. Therefore, 
the ability of a particular class of wave function to accurately 
reproduce the electric dipole moments of the fragments may be 
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some indication of its ability to predict their interaction 
energy. In Table IV are given theoretical and experimental 

values of the dipole moments of OH, H0^» and the HO HO^ 

complex. The calculated moments for OH and HO^ are too large by 
approximately 9 and 3X , respectively. If the interaction between 
the radicals was purely dipolar, the binding energy would be 
overeat imated by about 1 2X . This serves as an indication that 
the interaction energy may be slightly overestimated, but should 
not be greatly in error. 

It has been our observation that well-correlated wave 
functions frequently have less charge separation and smaller 
dipole moments than their uncorrelated < SCF ) counterparts . This 
observation is consistent with the results for two of the species 
in Table IV, but is the opposite of that for FO^. The HO^, wave 
function indicates a considerable charge transfer from the 
hydrogen to the central oxygen atom, as i3 shown in the table by 

44 

the Mulliken populations . In the Cl wave function the charge 
has been somewhat more evenly distributed over the two oxygen 
atoms than at the SCF level. Contributing to this redistr ibut ion 
is an important configuration (c=0.088) that is generated from 
con f 1 gurat i on (5) by the single excitation la” -> 2a”. This 
excitation moves an electron from a '•-type orbital that is mostly 
on the central oxygen to one that is primarily on the terminal 
oxygen atom. The use of a well-correlated wave function has 
resulted in a more even charge distribution between the oxygen 
atoms in HO^, but it is accomplished in a way that moves electron 
density to the terminal oxygen atom, increasing the overall 
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charge separation and dipole moment. The SCF dipole moment of 
H0 q actually agrees rather well with the experimental value. 

That situation appears to be the result of a fortuitous 
cancellation between the error due to the U3e of a limited < DZP ) 
basis set and thet due to use of an uncorrelated (SCF) wave 
f unction . 

Comparison of the electron populations for the complex and 
isolated radicals supports our earlier observation that the 
radicals change little upon binding. Table IV shows a net 
transfer of only 0.03 electrons from OH to HO^ upon complexat ion . 
A one to one comparison was made of the heavily occupied natural 
orbitals belonging to the isolated and complexed species. All 

except two of the HO HO^ natural orbitals have a very close 

resemblance to either an OH or HO^ orbital. However, two of the 
NO’s, both with populations of about 1.98, showed considerable 
change. In Fig. 2 are shown contour plots of the NOs that 
primarily represent the bond of HO^ for the isolated radical 

and the hydrogen bonded complex. It is easily seen from the 
figure that the negative end of the hydroxyl dipole has caused 
some of the charge density to push from the H-0 bonding region 
back onto the terminal oxygen atom. The population analysis for 
this orbital shows a transfer of about 0.4 electron from hydrogen 
and the central oxygen in H0 2 to the terminal oxvgen. The nodal 
structure of this orbital is such that this introduces some 0-0 
TT-type antibonding character. A quite different situation is 
found if the NO corresponding to the p-type lone pair on the 
terminal oxygen atom of HO^ is examined (Fig. 2). 


In this case 
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complexation with the OH radical causes this electron pair to 
delocalize significantly, transferring approx imately 0.35 
electron onto the central oxygen and hydrogen atoms. The nodal 
pattern of this orbital causes this transfer to introduce T-type 
0-0 bonding character. Thus, while individual natural orbitals 
undergo significant changes upon formation of the hydrogen bonded 
complex, they compensate in such a way that the total description 
changes only a little. It should be noted that, since this pair 
of natural orbitals is nearly degenerate with respect to 
occupation numbers, a more complete Cl description could result 
in their further mixing. 

The electric dipole moment calculated for the HO H0^ 

complex at the Cl level is 2.17 Debye. It is oriented to point 
roughly from the central oxygen atom in HO^ (negative end) toward 
the midpoint of the hydroxyl bond (positive end). This vector 
makes an angle of 72.1° with the 0-0 bond in H0,p. 

4. Comparison with H^O and H0^ dimers. 

The hydrogen bond in the water dimer has received a great 
deal of experimental and theoretical attention in recent years 
and is quite well characterized . Although this bond is weak when 
compared to the entire range of hydrogen bonds, it is an 
appropriate subject for comparison to the present system because 
both species involve bonding with the OH noiety. In Table V are 

54 5 i 55 

shown experimental and theoretical" ’ values of the water 
dimer electronic binding energy as well as those that we have 
calculated for the HO HO^ system. These latter values include 


the BSSE correction described above. 


The data in Table V suggest 
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that the hydrogen bond in HO HO^ is about 807 . as strong as that 

in the water dimer. The numbers reported for HO HO^ might 

decrease a few tenths of a kcal/mol if the DZP basis set were to 
be expanded considerably. This expectation is based on the 
observation that the present basis set yields dipole moments for 
the fragments that are a few percent too large. 

The experimental value^^ of the oxygen-oxygen distance in 

is 2.98 A and the geometry of the hydrogen bond ( 0-H 0) 

is linear to within 1-2 degrees. The equilibrium HO distance in 
water (0.9572 Al was used for the monomer in the above 

analysis^, yielding an experimental H 0 distance of 2.02 A. 

In the present case, the H 0 distance is only 0.02 A shorter 

than this value, although the 0-0 distance (2.85 A) is 0.13 A 
shorter than in the water dimer. This decrease is mostly the 

result of the hydrogen bond being bent in HO HO^, rather than 

the»-e being a large decrease in the H 0 or 0-H distances. The 

principal difference between the hydrogen bond in the present 
case and that in the water dimer is the significant deviation 
from linearity. It is noted, however, that the difference in SCF 
energy between the minimum energy structure (Fig. 1 1 and one 

formed by straightening the 0-H 0 linkage is only about 0.5 

kcal/mol . 

It is interesting to compare the hydrogen bonded HO HO^ 

system with the HO^ dimer that has recently been studied 
57 

theoretically. The HO^ dimer is considered a possible 
intermediate in the self recombinat ion of H0 o radicals, a 
reaction that also shows complex kinetic behavior. The 
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theoretical study , carried out at the SCF-DZP level, found a 
minimum energy conformation that was planar and cyclic with two 
equivalent hydrogen bonds. The total binding energy of the dimer 
was predicted to be about 4.9 kcal/mol. Although this is about 
the same stabilization energy as in the water dimer, there are 
two hydrogen bonds formed, and individually they must be 
considered very weak. Consistent with this interpretation, the 

H 0 bond lengths were found to be 2.19 A, or over 8% longer 

than in the water dimer. 

When compared with the HO^ dimer structure, that of HO HO^ 

is seen to be quite different. Although the general shape of the 

HO HO^ complex is suggestive of a quasi-cyclic nature, the 

H 0 distances do not support such an interpretation . The 

shorter H 0 distance is nearly the same as found in the water 

dimer, and the other one that would correspond to a second 

hydrogen bond is much too long, at 2.77 . Also, the HO HO^ 

potential surface is quite flat and very little energy (~0.5 
kcal/mol) is required to distort the complex away from the 
"cyclic" form. This is easily provided by thermal collisions at 
room temperature. 

Although the experimentalist is most likely to measure the 
room temperature enthalpy change of Reaction (7), it is not 
possible to estimate it quantitatively because some of the 
vibrational frequencies of the complex are unknown. This becomes 
especially critical when computing the room temperature energy of 
the low frequency inter-radical vibrations and hindered 


rotations . 


However, because the hydrogen-bond is similar to that 
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in the water dimer, it may be useful to treat these frequencies 

54 

as approximately equal in the two species. Both Curtiss et al . 

55 

and DelBene et al . have recently converted the enthalpy change 

for water dimer formation to the electronic energy change in 

order to compare experimental and theoretical results. In this 

context, the main differences between water dimer formation and 

Reaction (7) are that in the present case two degrees of 

rotational freedom are lost instead of three and there is one 

less low frequency vibrational mode gained. At room temperature 

these two effects are largely offsetting, allowing us to use the 

published water dimer result as a rough estimate. A correction 

55 

of ♦1.7 kcal/mol was found for the water dimer at 298 K. 

Applying this correction to the C I ( SDQ ) data in Table V yields an 
estimated enthalpy change for Reaction (7) at 298 K of -2.9 
kcal/mol . 

5. Other electronic states and isomers 

Calculations were carried out at the conformation of Fig. 1 

for three other electronic states of the complex. A ^A’ state 

can be made from configuration (8) by recoupling the two unpaired 

spins to yield a singlet state. The unpaired electrons are 

heavily localized in p -type orbitals that are centered 

approximately 3 A apart, and their coupling energy is expected to 

be small. In fact, the 1 A ' state was found to be 0.2 mh below 
3 

the A' state at the C I ( SDQ > level. This small energy difference 
is beyond the precis* "'n of the computational procedure, and the 
two states should be considered degenerate. There is also a pair 
of 1 A” and ^A” states that can be represented by the excitation 
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11a* > 2a” relative to configuration (8). Physically, these 

states correspond to placing the unpaired electron of the OH 

radical in an in-plane molecular orbital. The 1 A" and ^A” states 

3 

were found to lie 6-8 mh above the A’ ground state. It is 
emphasized that these energy gaps were calculated using the "A' 
minimum energy geometry, and no attempt has been made to explore 
the potential energy surfaces of the excited states. 

SCF calculations were performed for the ^A * state at 
conformations that would be reasonable choices for hydrogen 
bonded complexes with the hydroxyl radical acting as proton donor 
to either the central or terminal oxygen atom of the hydroperoxyl 
radical. Such conformations had already been rejected as energy 
minima by the gradient search procedure, and the energies at 
these conformations were found to be 4-5 mh above that for the 
structure in Fig. 1 . No attempt was made to determine the 
correlation energy contribution to these numbers, but a large 
differential effect would not be expected. Our conclusion is 
that the other conformations are probably not thermally 
accessible to a complex that has been stabilized at the low 
energy structure. It is clear, however, that there are large 
parts of the potential surface at relatively low energies and 
that, until stabilized, a hydrogen bonded complex will have a 
relatively large amount of phase space available to it. 
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IV. CONCLUSIONS 

This study has investigated possible intermediates in the 

reaction between hydroxyl and hydroperoxyl radicals. Two forms 

of the proposed intermediate were investigated. A covalently 

bonded form with the structure HOOOH is calculated to have an 

electronic binding energy of approximatel y 26 kcal/mol. 

Reasonable correction for basis set incompleteness would increase 

that energy by 1-2 kcal/mol. This value yields a predicted room 

temperature enthalpy change for Reaction (2) of -24.8 kcal/mol, 

which is aDproximately 5 kcal/mol less negative than a prediction 

made using group additivity considerations to estimate the heat 
25 

of formation of the complex. In addition, this study has 

identified a hydrogen bonded complex that has a binding energy of 

about 4.7 kcal/mol. There is no direct experimental or 

theoretical data for comparison, but we note that the bonding is 

similar to that found in the water dimer and considerably 

57 

different than predicted" for the HO^ dimer. 

The covalently bonded intermediate has symmetry and is in 

1 

a A electronic state, which cannot be correlated to ground state 

1 3 - 

products [H_0 ( A. ) and 0_( £ ) D via a spin allowed mechanism, 

c 1 2 g 

However, the initial product states in Reaction (1 1 are unknown. 

1 1 4 

Both of the excited-state oxygen channels ( A and £ 1 are 

g g 

energetically accessible, and cannot be ruled out as the primary 
reaction paths. The HOOOH 1 * 3 A state could be correlated with 
these asymptotes. 

The lowest energy form of the hydrogen-bonded intermediate 

3 1 

is found to be planar with A* and A' states that are 
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essentially degenerate. If the reaction were constrained to 

proceed along a planar path, these states would have incorrect 

symmetry for correlation with the ground state products ( 3 A" ) . 

The state could be correlated with the low-lying singlet 

excited states of the products, but the ' state must be 

correlated with excited triplet channels that are not 

energetical ly accessible. There is no ^ priori reason, however, 

for requiring the reaction to proceed with symmetry. If, 

during the course of the reactive event the molecule distorted 

away from planarity, the ground state triplet could mix with 

3 

other states of the correct symmetry ( A** ) to yield ground state 

products. Although geometry optimization was only carried out 
for the ^A ’ state, these calculations indicate that there is a 
low lying set of "* A’* and ^A" states that are very slightly bound 
with respect to the reactant radicals. 

As stated in the introduction, our results are not 

23 

consistent with the theoretical results of Rao et al. These 
authors have calculated the stability of various H^O^ forms 
relative to OH ♦ HO^ using smaller basis sets (ST0-3G, 4-31C) at 
the SCF level. They conclude that the hydrogen bonded form is 
not stable and that the covalent HOOOH intermediate has a 
stability of 6.3 kcal/mol . While our best Cl calculations 
clearly show HOOOH to be stable, tho SCF results would predict 
that it is unbound by about 10 kca/mol. To explore this 
disagreement further, we repeated the SCF calculations using the 
basis set (4-31G) and conformations given in Ref. 23. At this 
level our results would have predicted HOOOH to be even more 


unstable than with the larger basis set. The energies we 

obtained for the OH and HO.-, fragments were about 3 mh higher than 

reported in Ref. 23. Since we employed RHF calculations , this 

small discrepancy would be expected if the previous results were 

obtained with a UHF procedure. For HOOOH, however, our 4-31C 

energy is -225.2123 hartree, about 48 mh above that of Ref. 23. 

There is no obvious explanation of this large discrepancy. It 

20 22 24 

should be noted that several previously published reports ’ 

include SCFC4-31G) treatments of HOOOH, and that all of these 

23 

results are in disagreement with that of Rao et al . 

It is equally difficult to understand the prediction of Ref. 

23 that the hydrogen bonded form of H^O^ would not be stable. It 

is common for unpolarized basis sets such as the 4-31G to give 

dipole moments that are too large for the unbound fragments. In 

a hydrogen-bonding situation, this defect usually leads to 

binding energies that are much too large. For example, the 4-31G 

5 8 

set has been shown to yield an SCF stabilization energy for the 
water dimer of 8.2 kcal/mol without correction for BSSE . Some of 

our preliminary SCF gradient calculations on HO HO^ were 

carried out with a 4-31G basis set. While the geometry was never 
completely optimized, energies were obtained that were lower than 
the fragment energies published in Ref. 23 and would have 
indicated binding. The only simple explanation for this 
discrepancy is that tne geometry search procedures of Ref. 23 did 
not sample regions of the potential surface sufficiently close to 


the energy minimum. 
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One other possibility that conies to mind is that the authors 
of Ref. 23 have simply misidenti f ied their calculated values for 
the covalent and hydrogen bonded forms. The energy they have 
published for the covalent form (-225.260 au ) would not be an 
unreasonable 4-31G UHF value for the hydrogen bonded complex. If 
this change in identification were made, their results would be 
in qualitative agreement with the present ones. 

Of the two intermediates studied here, the hydrogen bonded 
form is not as stable as is usually suggested in proposed 
mechanisms 1 ^ for Reaction (1 ). However, we note that the 
hydrogen bonded part of the potential surface is rather flat and 
that a large amount of phase space would appear to be accessible 
to the reactant radicals during complex formation. In addition 
to the nearly degenerate singlet and triplet A* states, the low 
lying singlet and triplet A” states have energies below that of 
the reactants and can be expected to contribute to the dynamics. 

The stability of the covalently bonded intermediate is in the 

1 

usual range. The test of whether or not either of these species 
can serve as intermediates in a mechanism that explains the 
experimental data is in the application of kinetics models to the 
problem. In order to do this, we must know the barriers to 
product formation for each of these intermediates as well as the 
vibrational frequencies in the complexes and their transition 
states. Work in both of these directions is proceeding in these 
1 abora tor i es . 
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TABLE I . 
molecule 


Total and binding energies (in hartrees ) of the HOOOH 



SC F 

CI( SB ) 

C 

o 

C I ( SPG ) 

E( OH ) 

-75.4063 

-75.5630 

0.9775 

-75.5700 

EC H0 2 ) 

-150.2121 

-150.5377 

0.9525 

-150.5678 

E( OH ) ♦ E(H0 2 ) a 

-225.6184 

-226.1 007 

. . . 

-226.1 378 

E( OH ♦ H0 2 > b 

-225.61 84 

-226.0755 

0.9424 

-226.1 267 

E( HOOOH ) 

-225.601 9 

-226.1 079 

0.9380 

-226. 1 688 

c 

Binding Energy 

-0.0165 

♦0 . 0324 

. . . 

♦0 . 0421 


a. Calculated separately for each radical. 

b. Calculated as supermolecule with radicals separated by 500 a . 


c . 


Calculated 


as 


Line 5 


Line 4 
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Table II. Reaction energies for formation of HOOOH . 




Energy(kcal/mol ) 

Description 3 


AE 

e 

-26.4 

Electronic energies only 

AE 

o 

-22.2 

''ibrational zero pt . energy 

AE ,b 

o 

-23.5 

ad hoc correction for basis 
set incompleteness. 

AE ’ 

298 

-24.3 

Vib, rot, and trans energy 
at 298 K . 

AH * 

298 

-24.8 

APV work term < -RT ) . 


a. Indicates what contribution has been included to obtain entry 
from previous one. 

b. The prime indicates that the basis set correction is included. 
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TABLE III- Total and binding energies (in hartrees ) for the 
HO HO^ complex. 9 



SCF 

Cl ( SD ) 

C 

o 

CICSDQ ) 

E( OH ) 

-75.4065 

-75.5627 

0 . 9778 

-75.5696 

EC H0 2 > 

-150.21 33 

-150.5365 

0.9538 

-150.5657 

EC OH) ♦ E(H0 2 > b 

-225.61 98 

-226.0992 

. . . 

-226.1 353 

EC OH ♦ H0 2 ) C 

-225.61 98 

-226.0742 

0.9435 

-226.1 240 

EC HO H0 q ) 

-225.6272 

-226.0842 

0.9457 

-226.1 325 

Binding Energy^ 

0.0074 

0.01 00 

• ■ ■ 

0 .0085 


a . 

b. 


c . 
d . 


Calculated 

Calculated 

Calculated 

Calculated 


using UHF gradient conformations, 
separately for each radical, 
as supermolecule with radical separated 
as Line 5 - Line 4. 


by 50 0 
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TABLE IV. Electric dipole moment and Mulliken populations for 
the HO, HO.,, and HO HO^ species. 

a 

Dipole Moment Populations 



( Debye ) 

HI 

01 

02 

H2 

03 

'2 

SCF 

2.08 

0.64 

8.30 

8.06 

• ■ • 

• • • 

Cl 

2.16 

0.65 

8.24 

8.10 

• • • 

■ • • 

EXP 

2. 09 b 

■ • • 

• ■ • 

• • • 

. . . 

. . . 

SCF 

1 .91 

• • • 

9 m m 

• • • 

0.66 

8.34 

Cl 

1 .82 

• • • 

• • • 

• • • 

0.68 

8.32 

EXP 

1.67° 

• • • 

• • • 

. . . 



1 — ho 2 

SCF 

2.28 

0.^3 

8.29 

8.10 

0.63 

8.35 

Cl 

2.17 

0.64 

8.25 

8.14 

0.64 

8.33 


a. H0-, and OH are labelled as follows: HI-01-02 and H2-03. 

b. Reference 59. 

c . Ref erence 60 . 


m 


T 
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Table V. Electronic binding 
bonded species. 

energies for three 

hydrogen 

E( kcal/mol ) 

Method 

Reference 

(HO) 

2 2 

5.44 

Expt 1 

55 

3.6 

SCF 

51 

5.5 

Cl 

51 

4.3 

SCF 

56 

5.4 

MBPT 

56 

(H0 2 ) 2 

4.9 

SCF 

57 

HO H0 o 

4.0 

SCF 

This study 

4.7 

C I ( SDQ ) 

This study 


I 

» 


) 
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FIGURE CAPTIONS 

1. Minimum energy conformation of hydrogen bonded HO HO^ 

complex, determined at the UHF level. 

2. The H-0 sigma bonding natural orbital in isolated HO^ ( top "> 
and in the hydrogen bonded complex (bottom). Contour intervals 
are 10% of the function’s range. 

3. The in-plane lone pair natural orbital in isolated HO^ (top) 
and in the hydrogen bonded complex (bottom). Contour intervals 
are 10% of the function’s range. 
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a = 104 . 3 * 
ft = 147 . 5 * 
)'= 104 . 5 * 


\ 



2.003 


\ 


\ 

\ 


A 



A 





